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Capacitance  and  resistance  are  the  fundamental  electrical  parameters  used  to  evaluate  the  electrical 
characteristics  of  a  supercapacitor,  namely  the  dynamic  voltage  response,  energy  capacity,  state  of  charge 
and  health  condition.  In  the  British  Standards  EN62391  and  EN62576,  the  constant  capacitance  method 
can  be  further  improved  with  a  differential  capacitance  that  more  accurately  describes  the  dynamic 
voltage  response  of  supercapacitors.  This  paper  presents  a  novel  bivariate  quadratic  based  method  to 
model  the  dynamic  voltage  response  of  supercapacitors  under  high  current  charge-discharge  cycling, 
and  to  enable  the  derivation  of  the  differential  capacitance  and  energy  capacity  directly  from  terminal 
measurements,  i.e.  voltage  and  current,  rather  than  from  multiple  pulsed-current  or  excitation  signal 
tests  across  different  bias  levels.  The  estimation  results  the  author  achieves  are  in  close  agreement  with 
experimental  measurements,  within  a  relative  error  of  0.2%,  at  various  high  current  levels  (25—200  A), 
more  accurate  than  the  constant  capacitance  method  (4-7%).  The  archival  value  of  this  paper  is  the 
introduction  of  an  improved  quantification  method  for  the  electrical  characteristics  of  supercapacitors, 
and  the  disclosure  of  the  distinct  properties  of  supercapacitors:  the  nonlinear  capacitance- voltage 
characteristic,  capacitance  variation  between  charging  and  discharging,  and  distribution  of  energy  ca¬ 
pacity  across  the  operating  voltage  window. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  the  electrical  modelling  of  supercapacitors,  capacitance  and 
resistance  are  the  fundamental  electrical  parameters  used  to  eval¬ 
uate  the  electrical  characteristics  of  a  supercapacitor,  namely  the 
dynamic  voltage  response,  energy  capacity,  energy  loss,  state  of 
charge  and  health  condition.  The  capacitance  of  supercapacitors 
(electrochemical  capacitors)  consists  of  two  main  charge  storage 
principles:  double-layer  capacitance  and  pseudocapacitance  [1,2], 
depending  on  the  electrode  material  used.  In  double-layer 
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capacitance,  charge  is  stored  electrostatically  in  the  interface  be¬ 
tween  an  electrode  and  an  electrolyte.  When  an  external  voltage  is 
applied  on  one  electrode,  the  counter-ions  of  the  electrolyte  are 
attracted  to  the  surface  of  the  electrode,  forming  two  layers  of 
charge  that  separated  by  a  molecular  dielectric  layer  (mono-layer  of 
solvent  molecules)  at  the  electrode-electrolyte  interface.  In  pseu¬ 
docapacitance,  the  storage  mechanism  is  a  reversible  Faradaic 
charge-transfer  process,  achieved  by  the  adhesion  of  electrolyte 
ions  onto  the  surface  of  electrode  (made  up  of  metal  oxides  or 
conducting  polymers)  via  electrosorption,  intercalation  or  reduc¬ 
tion-oxidation  reactions. 

The  capacitance  of  a  supercapacitor  is  susceptible  to  changes  in 
operating  conditions  such  as  voltage,  frequency,  temperature,  and 
charge  propagation  [3-5  ,  mainly  related  to  the  mobility  of  charge 
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during  operation.  The  voltage  dependency  of  the  device  capaci¬ 
tance  can  be  explained  by  the  effect  of  the  applied  electric  field  on 
the  permittivity  of  the  molecular  dielectric,  the  separation  of  the 
charge  layer  and  the  molarity  of  electrolyte.  Besides,  the  capaci¬ 
tance  is  proportional  to  the  effective  area  of  the  charge  layer  that 
formed  at  the  electrode-electrolyte  interface  during  operation,  not 
the  total  surface  area  of  typically  highly  porous  electrode.  This 
effective  area  depends  on  the  accessibility  of  electrode  pores  and 
the  propagation  of  electrolyte  ions  into  the  different  pore  depths/ 
sizes  of  electrode.  The  propagation  of  ions  takes  longer  time  to 
reach  long  narrow  pores  than  short  wide  pores,  hence  resulting  in 
higher  capacitance  (due  to  larger  effective  area)  at  lower  frequency 
operation.  This  explains  both  the  charge  and  frequency  effects  on 
the  capacitance. 

In  the  British  Standards  EN62391  [6]  and  EN62576  [7]  (iden¬ 
tical  to  IEC62391  and  IEC62576),  the  capacitance  of  super¬ 
capacitors  is  quantified  as  a  constant  capacitance  by  using  the 
constant  current  discharge  test  method,  without  considering  the 
voltage  and  charge  dependencies  of  the  supercapacitor.  This 
constant  capacitance  based  identification  method  is  found  to  be 
inadequate  to  precisely  represent  the  dynamic  terminal  behaviour 
of  supercapacitors  [8-10]. 

The  term  differential  capacitance  is  the  measurement  of 
capacitance  as  a  function  of  voltage.  Differential  capacitance  is 
commonly  introduced  within  the  electrical  models  of  super¬ 
capacitor  to  describe  the  voltage  dependency  of  capacitance.  In  the 
deployment  of  differential  capacitance  in  various  resistor-capac¬ 
itor  equivalent  circuit  models,  the  main  challenge  is  obviously  to 
accurately  establish  the  relationship  of  capacitance  versus  voltage, 
which  can  be  identified  from  either  the  frequency  response  to  a 
small  excitation  signal  (typically  currents  of  a  few  mA),  e.g.  in  tests 
using  electrochemical  impedance  spectroscopy  (EIS)  [4,11],  or  the 
voltage  response  to  a  controlled  pulse-current  (in  tens  of  Amperes) 
in  studies  using  the  constant  current  test  method  [12-14  repeat¬ 
edly  across  different  bias  levels.  The  differential  capacitance  can 
then  be  modelled  as  a  linear  function  [4,12,14-16  ,  a  quadratic 
function  [11  ,  or  even  a  quartic  function  [17  .  The  order  of  function 
relates  to  the  size  of  the  experimental  dataset  (number  of  bias 
levels)  used  in  defining  the  relationship  of  capacitance  versus 
voltage.  Studies  with  small  datasets  consume  less  experimental 
time  and  generally  allow  simple  analysis  compared  to  studies  with 
large  datasets;  however  the  precision  of  modelling  could  be 
comprised. 

Apart  from  the  linear  capacitance-voltage  based  equivalent 
circuit  models,  several  analytical  models  have  been  proposed  to 
model  the  dynamic  behaviour  of  supercapacitors.  These  models  are 
ladder  network  model  with  a  complex  pore  impedance  block  18], 
transmission  line  model  with  the  hybridisation  of  temporal  and 
frequency  approaches  11],  artificial  neural  network  based  models 
[19—21  ,  and  fractional  non-linear  model  [22  .  In  these  models,  the 
parameters  identification  often  involves  intensive  experimentation 
and  computation  to  achieve  good  approximation  results. 

This  paper  presents  a  highly  accurate  and  practical  method  to 
model  the  dynamic  voltage  response,  and  quantify  the  differential 
capacitance  and  energy  capacity  of  supercapacitors  under  high 
current  operations  (up  to  200  A)  that  typically  experienced  in  real- 
life  applications.  A  bivariate  quadratic  function  is  developed  and 
applied  for  the  first  time  to  represent  the  dynamic  charge-voltage 
characteristic  of  supercapacitors.  This  proposed  method  allows  the 
differential  capacitance  and  energy  capacity  of  supercapacitors  to 
be  derived  directly  from  the  terminal  measurements,  i.e.  voltage 
and  current,  rather  than  from  multiple  pulse-current  or  excitation 
signal  tests  at  different  bias  levels  as  discussed  previously.  Hence, 
the  bivariate  quadratic  based  method  can  eliminate  the  trade-off 
between  accuracy  and  practicality  (i.e.  simplicity  of  testing)  in 


identifying  the  capacitance-voltage  characteristic  for  the  electrical 
modelling  of  supercapacitors. 

2.  Experimental 

This  study  aims  to  understand  the  dynamic  behaviour  of 
supercapacitors  during  the  high  power/current  operations  often 
experienced  in  real-life  applications.  In  the  laboratory,  these 
operating  conditions  are  recreated  as  constant  current  charging¬ 
discharging  experiments  by  using  a  programmable  high  current 
test  system.  The  high  current  test  system  10]  developed  by  the 
authors,  supports  three  types  of  charge/discharge/cycling  current 
profile  -  constant  current,  sinusoidal  current,  and  customised 
current  cycles,  at  current  levels  up  to  250  A.  A  data  acquisition 
device,  based  on  a  National  Instrument  USB-6218,  is  deployed  to 
record  the  voltage,  current  and  temperature  measurements  of  the 
supercapacitor  under  test. 

High  current  charging-discharging  experiments  were  con¬ 
ducted  on  Maxwell’s  Boostcap  supercapacitors.  The  specifications 
of  these  double-layer  based  supercapacitors  are  an  equivalent  se¬ 
ries  resistance  (ESR)  of  0.8  mQ,  a  nominal  capacitance  of  650  F,  and 
a  rated  voltage  of  2.7  V.  From  a  constant  current  charging-dis¬ 
charging  experiment  of  200  A,  the  voltage  response  of  a  650  F 
supercapacitor  cell  operated  between  0  and  2.7  V  is  as  depicted  in 
Fig.  1.  The  current  signal  is  the  effective  square  wave,  with  a  cor¬ 
responding  triangular  voltage  response.  The  ‘charge’  currents  are 
represented  by  negative  values  whilst  the  ‘discharge’  currents  are 
represented  by  positive  values.  Both  the  voltage  and  current  time 
series  data  are  subsequently  used  as  the  input  data  for  the  proposed 
numerical  method,  described  in  following  section. 

3.  Bivariate  quadratic  method 

3.1.  Simple  resistor-capacitor  circuit 

For  practical  purposes,  a  simple  series  resistor-capacitor  (RC) 
equivalent  circuit  is  applied  to  represent  the  dynamic  voltage 
response  of  supercapacitors  during  fast  charging-discharging.  The 
circuit  parameters  are  stated  as  the  equivalent  series  resistance 
(ESR)  and  differential  capacitance,  as  shown  in  Fig.  2.  In  terms  of 
circuit  notation,  the  terminal  voltage  of  the  supercapacitor  is 
denoted  as  Use,  the  current  of  the  supercapacitor  as  /sc,  the  voltage 
across  the  ESR  as  Vr,  and  the  differential  capacitor  voltage  is  Vo 


Fig.  1.  The  voltage  response  of  a  supercapacitor  under  charge-discharge  cycling  at 
200  A. 
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Fig.  2.  Simple  RC  equivalent  circuit  of  a  supercapacitor. 

Within  the  context  of  this  paper,  the  “capacitor”  term  refers  to  the 
differential  capacitance,  not  the  entire  supercapacitor  cell. 


3.2.2.  Equivalent  series  resistance 

The  voltage-time  characteristic  of  a  supercapacitor  during  a 
constant  current  charging-discharging  cycle  is  illustrated  in  Fig.  3. 
The  ESR  parameter  can  be  identified  from  the  abrupt  voltage 
change  (rise  or  fall)  at  the  start  or  end  of  the  charging  or  discharging 
state.  This  abrupt  voltage  change  is  due  to  the  change  in  voltage 
across  the  ESR,  Vr,  when  there  is  a  step  change  in  current,  A/so 

ESR  =  ff-  (1) 

msc 

From  Figs.  1  and  3,  the  abrupt  voltage  change  at  the  start  of  the 
discharging  state  (time  t2)  is  twice  the  level  at  the  start  of  charging 
(time  ti)  due  to  a  twofold  increase  in  the  step  change  in  current 
(A/sc  =  2  x  /sc).  Thus,  the  ESR  determined  at  times  t\  and  t2  should 
be  close  in  value  to  each  other.  Vm[n  and  l/max  are  the  minimum  and 
maximum  voltage  points  defined  during  the  constant  current 
cycling  experiment. 


3 A. 2.  Differential  capacitance 

A  linear  capacitor,  by  definition,  exhibits  a  linear  increase  in 
voltage,  V ,  with  respect  to  charge,  Q,  implying  its  capacitance,  C,  is 
always  constant. 


Q  =  dQ 

V  dV 


A  supercapacitor  (a  form  of  non-linear  capacitor)  exhibits  a 
differential  increase  in  voltage,  V \  with  respect  to  charge,  Q.  The 
differential  capacitance,  Cdiff,  is  obtained  from: 


By  referring  to  (3),  the  differential  capacitance  of  a  super¬ 
capacitor  can  be  determined  from  the  derivative  of  its  charge- 
voltage  characteristic.  In  order  to  establish  the  relationship  of 
voltage  versus  charge,  the  charge  exchanged  during  the  experiment 
is  calculated  by  using  the  measured  current  time  series  data.  At  a 
constant  current  /sc,  the  charge  exchanged  during  charging  (be¬ 
tween  t\  and  t2),  Qchg,  and  charge  exchanged  during  discharging 
(between  t2  and  £3),  Qdchg,  are  determined  from: 

h 

Qchg  =  J  hcdt  (4) 


3.2.  Bivariate  quadratic  equation 

By  considering  the  voltage  and  charge  dependencies  of  the 
capacitance,  the  authors  propose  a  new  bivariate  quadratic  based 
method  to  approximate  the  charge-voltage  characteristic  of  a 
supercapacitor.  The  derivative  and  integral  of  the  proposed 
approximation  are  used  to  quantify  the  differential  capacitance  and 
energy  capacity  properties  of  the  supercapacitor  respectively.  The 
bivariate  quadratic  equation  is  therefore  of  the  following  form: 

Ax2  +  Bxy  +  Cy 2  +  Dx  +  Ey  =  1  (6) 

where  x  relates  to  charge,  Q,  y  relates  to  voltage,  V,  and  the  co¬ 
efficients,  A  to  E,  require  to  be  determined. 

The  proposed  approximation,  developed  with  the  Matlab®  nu¬ 
merical  computing  environment,  comprises  four  simple  steps  to 
accurately  represent  the  charge-voltage  characteristic  of  a  super¬ 
capacitor.  The  procedures  are  as  follows: 


3.2 A.  Step  1:  select  five  data  points 

With  five  unknown  coefficients  in  Equation  (6),  five  data  points, 
Ptl  to  Pt5,  are  required  to  be  chosen  from  the  charge-voltage 
characteristic  so  that  these  points  are  equally  spaced  along  the 
charge,  Q,  axis,  as  shown  in  Fig.  4. 


Fig.  3.  The  voltage-time  characteristic  of  a  supercapacitor  under  constant  current 
cycling. 


Fig.  4.  The  charge-voltage  characteristic  of  supercapacitor  under  constant  current 
cycling. 
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3.2.2.  Step  2:  translate  coordinate  plane 

Numerical  computation  on  the  five  unknown  coefficients  of  A  to 
E  could  result  in  unlimited  possible  values,  hence  in  order  to 
simplify  the  computation  process,  the  charge-voltage  plane  (Q— V 
plane)  is  translated  into  the  x-y  plane  so  that  the  two  data  points  of 
Ptl  and  Pt5  intersect  at  the  new  horizontal,  x,  and  vertical,  y,  axes 
respectively. 

During  the  charging  state,  the  translating  equations  are: 

X  =  o-  Qchg  (7) 

y  =  V  —  Vchgl  (8) 

During  the  discharging  state,  the  translating  equations  are: 

x  =  Q  (9) 

y  =  (10) 


3.2.3.  Step  3:  determine  the  coefficients  (A  to  E)  of  the  bivariate 
quadratic  equation 

By  substituting  the  translated  data  points  (in  x-y  plane)  into 
Equation  (6),  the  resulting  equations  can  be  written  in  the  matrix 
form  as  follow: 


'*i  *iyi  y\  yr 

"A" 

"1" 

x\  x2y2  y\  x2  y2 

B 

1 

*3  *3y3  *3  y3 

• 

C 

— 

1 

*4  *4y4  y\  *4  y4 

D 

1 

x\  x5y5  y\  x5  y5_ 

E 

1 

The  coefficients  of  A  to  E  can  be  determined  by  using  the  matrix 
inversion  method. 


A 

B 

C 

D 

E 


x\y\ 

y\ 

^1 

yr 

X2V2 

y\ 

X2 

yz 

y\ 

X3 

y3 

x4y4 

y\ 

X4 

y4 

x5y5 

y§ 

X5 

ys_ 

-1 

"i" 

1 

. 

1 

1 

1 

(12) 


3.2.4.  Step  4:  approximate  the  charge -voltage  characteristic  of 
super  capacitors 

For  a  quadratic  function/(y)  =  ay2  +  hy  +  c,  its  roots  (values  ofy 
which  fly)  =  0)  are 


-b±VK 
y  =  ^uT~ 


(13) 


where  the  discriminant,  A  =  b2  -  4 ac. 

Once  the  coefficients  of  A  to  E  are  identified,  the  approximation 
of  the  voltage  response  of  the  supercapacitor  with  respect  to  charge 
is  computed  according  to  the  root-finding  method  described  in 
Equation  (13). 

Re-arranging  Equation  (6)  given: 


y 2  + 


Bx  +  E\  Ax2  +  Dx  -  1  _ 

y  + - = - =  o 


c 


c 


The  roots  of  Equation  (14)  are: 


(14) 


\  (  [Bx  +  £\  /  Tbx  +  £\2  ./Ax2+Dx-1\\ 

y=2[-{—)±n—)  — c — )) 

(15) 

By  substituting  the  translating  Equations  (7)— (11)  into  Equation 
(15),  the  charge-voltage  characteristic  of  the  supercapacitor  is 
approximated  by  the  following  equations. 

During  the  charging  state: 


^sc(Q)  =  5 


Ax2  +  Dx  -  1 
C 


+  v< 


chgl 


(16) 


where  x  =  Q  -  Qchg. 

During  the  discharging  state: 


^sc(Q-) 


4(«!± DQ-1  ,  >  + 


(17) 


Next,  the  derivative  of  x  with  respect  toy  (dx/dy)  of  the  bivariate 
quadratic  equation  is  determined  via  implicit  differentiation. 


dx 


Bx  +  2Cy  +  E 
dy  2  Ax  +  By  +  D 


(18) 


As  the  translation  from  the  Q-V  coordinate  plane  to  the  x-y 
coordinate  plane  has  no  effect  on  the  gradient  of  the  polynomial, 
the  differential  capacitance  Cdiff  in  Equation  (3)  can  be  written 
as: 


Chi 


dQ.  dx 


diff 


dV'  dy 
During  the  charging  state: 


(19) 


Qliff_chg 


fi(Q-Qchg)+2  c(v-ychg,)+E 
2A(Q_-Qchg)+B(V-VchgA+D 


(20) 


During  the  discharging  state: 


Qliff_dchg 


BQ  +  2C(V  —  Umin)  +  E 
2  AQ  +  B(V-Vmin)+D 


(21) 


3.3.  Energy  capacity  of  a  supercapacitor 

At  a  constant  charge  current  /sc,  the  amount  of  energy  a 
supercapacitor  exchanged  (stored  in  charging  state  and  released  in 
discharging  state)  can  be  determined  by  calculating  the  area  under 
the  charge-voltage  curve,  as  shown  in  Fig.  5.  The  energy  dissipated 
in  the  ESR  is  denoted  as  Er,  the  energy  exchanged  in  the  differential 
capacitor  is  denoted  as  Eg  also  known  as  capacitor  energy,  and  the 
energy  remaining  inside  the  capacitor  before  the  charging  or  after 
the  discharging  is  denoted  as  E^e.  The  E^e  is  only  applicable  to  those 
experiments  with  a  non-zero  minimum  voltage  (or  end  voltage), 

^min* 
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Fig.  5.  The  energy  components  of  supercapacitor  under  constant  current  charging¬ 
discharging. 


The  quantification  of  capacitor  energy  is  achieved  by  applying 
the  trapezium  rule.  By  discretising  the  area  under  the  charge- 
voltage  curve  into  N  equally  spaced  panels,  the  capacitor  energy  is 
quantified  as  follows: 

During  the  charging  state,  the  dissipated  energy  in  ESR,  ER: 

Er  =  Vr_c hg  x  Qchg  (22) 

The  capacitor  energy,  Ec'. 


Ec. chg  - h  (l(Vsc(Qi )  +  Vsc(Qn+i  ))  +  E  VsdQk)  J  -  Vchgl  x  Qch 


g 

(23) 


where  spacing,  h  =  Qchg/Nand  Vsc(Q)  was  defined  in  Equation  (16). 

During  the  discharging  state,  the  remaining  energy  inside  the 
capacitor,  ERe: 

Er0  =  ^min  x  Qdchg  (24) 

The  capacitor  energy,  Ec : 


£c_dchg  ~  h  (  ^  (^Sc(Ql  )  +  ^Sc(Qn+1  ))  +  ^2^sc(0l<)  )  -  ^min  x  Qdchg 

(25) 


k= 2 


where  spacing,  h  =  Qdchg/Nand  Vsc(Q) was  defined  in  Equation  (17). 


3.3.1  Energy  equivalent  capacitance 

For  comparative  reasons  discussed  later,  the  following  energy 
equivalent  linear  capacitance  is  now  defined.  At  a  given  capacitor 
voltage  Vc,  the  energy  equivalent  capacitance,  Q,  of  a  linear 
capacitor  storing  the  same  amount  of  energy  as  the  differential 
capacitance  of  the  supercapacitor  is  given  as: 

2  Er 

=  (26) 


4.  Results  and  discussion 

In  this  section,  the  modelling  accuracy  of  the  proposed  bivariate 
quadratic  method  in  representing  the  dynamic  voltage  response  of 
supercapacitors  under  high  current  charging-discharging  is  eval¬ 
uated  and  compared  with  the  constant  capacitance  method.  The 
quantification  of  differential  capacitance  and  energy  capacity 


(capacitor  energy)  which  have  been  derived  directly  from  the  ter¬ 
minal  measurements,  are  presented  and  discussed. 


4 A.  Experimental  results 

The  electrical  characteristics  of  650  F  Maxwell  Boostcap  super¬ 
capacitors  were  evaluated  using  constant  current  charge- 
discharge  experiments  at  current  levels  of  25  A,  50  A,  100  A  and 
200  A.  In  these  experiments,  the  initially  uncharged  super¬ 
capacitors  were  cycled  between  a  minimum  voltage  of  0  V  and  a 
maximum  voltage  of  2.7  V.  At  the  end  of  each  cycle,  the  lowest 
achievable  voltage  level  was  limited  to  0.06  V  due  to  the  tolerance 
in  the  voltage  sensing  circuit  of  the  high  current  test  system.  The 
following  discussion  is  based  on  a  ‘fresh’,  i.e.  new,  650  F  Maxwell 
Boostcap  supercapacitor  (labelled  as  E22),  which  has  an  ESR  of 
0.8  mQ  and  a  nominal  capacitance  of  650  F. 

The  voltage  response  of  the  supercapacitor  E22  during  the  first 
cycle  of  the  constant  current  experiment  of  25  A  is  shown  in  Fig.  6. 
An  abrupt  voltage  change  of  0.012  V  responding  to  a  step  change  in 
current  of  24.9  A  is  recorded  at  the  start  of  charging  state  whilst  an 
abrupt  voltage  change  of  0.05  V  responding  to  a  step  change  in 
current  of  52.8  A  is  recorded  at  the  start  of  discharging  state.  The 
resultant  ESR  is  computed  to  be  0.47  mQ  in  the  charging  state 
whilst  0.95  mQ  in  the  discharging  state.  These  computed  ESR 
values  are  close  to  the  specified  ESR  of  0.8  mO.  During  the  charging 
state,  the  measurements  of  exchanged  charge,  internally  dissipated 
energy  and  capacitor  energy  are  1694.6  C,  20.0  J  and  2449.6  J 
respectively.  During  the  discharging  state,  the  measurements  of 
exchanged  charge,  remaining  energy  (due  to  non-zero  end  voltage) 
and  capacitor  energy  are  1617.6  C,  207.9  J  and  2124.8  J  respectively. 
Similar  measurements  were  obtained  from  constant  current  cycling 
experiments  carried  out  at  50  A  and  200  A,  and  were  then  compiled 
and  shown  in  Table  1. 

By  referring  to  the  experimental  results  in  Table  1,  the  ESR  was 
found  to  be  between  0.5  mQ  and  1.0  mQ,  close  to  the  specified 
value  of  0.8  mQ.  This  indicates  that  the  supercapacitor  E22  is  in 
pristine  condition.  Next  observation  is  that  an  increase  in  oper¬ 
ating  current  level  results  an  increase  in  the  dissipated  energy,  ER , 
and  a  decrease  in  the  exchanged  charge,  Q,  and  capacitor  energy, 
Ec,  of  the  supercapacitor.  This  can  be  attributed  to  the  increase  in 
the  voltage  drop  across  the  ESR  that  reduces  the  effective  working 
voltage  window  of  the  supercapacitor.  Next,  the  remaining  energy, 
ERe,  at  the  end  of  the  constant  current  cycle  is  caused  by  the  lowest 


Time,  t  (Seconds) 


to 

O.' 

<n 

Q- 

£ 

o 

to 


Fig.  6.  The  terminal  response  of  supercapacitor  E22  during  the  first  cycle  of  25  A 
cycling  experiment. 
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Table  1 

Experimental  results  (first  cycle)  of  supercapacitor  E22  at  different  current  levels. 


Current  /Sc  (A) 

Charging  state 

Discharging  state 

ESR  (mQ) 

Qchg  (C) 

£r(J) 

E-C c hg  (J) 

ESR  (mQ) 

Qdchg  (C) 

ERe  (J) 

Ec  dchg  (J) 

25 

0.47 

1694.6 

20.0 

2449.6 

0.95 

1617.6 

207.9 

2124.8 

50 

0.52 

1663.7 

43.2 

2368.0 

0.94 

1596.9 

95.8 

2102.1 

100 

0.68 

1639.6 

112.7 

2303.7 

1.03 

1532.3 

107.3 

1942.7 

200 

0.69 

1574.0 

220.4 

2175.4 

0.92 

1426.0 

92.7 

1695.1 

voltage  level  of  0.06  V  that  can  be  achieved  by  the  high  current 
test  system.  In  the  25  A  experiment,  the  FRe  is  twofold  the  values 
recorded  in  50-200  A  experiments  due  to  the  recorded  end 
voltage  of  0.129  V. 

4.2.  Charge -voltage  characteristic 

By  applying  the  procedures  described  in  Section.  3.2,  a  bivariate 
quadratic  equation  was  developed  from  five  experimental  data 
points  to  represent  the  charge-voltage  characteristic  of  super¬ 
capacitor  E22.  In  Figs.  7  and  8,  the  charge-voltage  characteristic 
during  25  A  charging  and  discharging  are  shown,  together  with  the 
approximation  results  based  on  the  constant  capacitance  method 
and  the  bivariate  quadratic  method.  The  experimental  data  is 
plotted  is  blue  (in  the  web  version);  the  constant  capacitance 
approximation  is  plotted  in  green;  and  the  proposed  bivariate 
quadratic  approximation  is  plotted  in  red.  In  both  figures,  the 
experimental  data  is  closely  represented  by  the  author’s  proposed 
model  compared  to  the  conventional  constant  capacitance  method. 
The  gradient  of  the  estimated  charge-voltage  curve  is  assessed  to 
quantify  the  differential  capacitance  of  supercapacitor  E22  in  Sec¬ 
tion.  4.3.  Another  observation  is  that  the  charging  phase  charge- 
voltage  curve  exhibits  larger  variation  in  the  gradient  than  the 
discharging  phase  charge-voltage  curve. 

4.3.  Differential  capacitance 

The  differential  capacitance  of  supercapacitor  E22  is  identified 
from  the  derivative  of  the  bivariate  quadratic  equation,  i.e.  Equa¬ 
tions  (20)  and  (21 ).  The  voltage  dependency  of  this  capacitance  at 
current  levels  of  25  A,  50  A,  100  A  and  200  A  is  shown  in  Fig.  9, 
with  the  capacitance  increasing  with  terminal  voltage.  The  rela¬ 
tionship  of  capacitance  versus  voltage  is  more  complex  than  a 
linear  function  which  is  suggested  in  some  characterisation 


Fig.  7.  The  charge-voltage  characteristic  of  supercapacitor  E22  under  25  A  charging. 


studies.  In  the  25  A  charging  state,  the  capacitance  is  found  to  be 
630.41  F  using  the  constant  capacitance  method  (Equation  (2)) 
whilst  in  the  range  of  456  F-749  F  using  the  author’s  proposed 
method.  In  the  25  A  discharging  state,  the  capacitance  is  641.53  F 
with  the  constant  capacitance  method  whilst  in  the  range  of 
557  F-696  F  with  the  proposed  method.  In  addition  to  the  voltage 
dependency  of  capacitance,  the  capacitance  variation  between 
charging  and  discharging  states  should  be  considered  in  electrical 
modelling  of  supercapacitors.  By  applying  the  differential  capac¬ 
itance  identified  from  the  proposed  method,  the  dynamic  voltage 
response  of  supercapacitor  E22  is  well  predicted,  as  shown  in 
Fig.  10. 

4.4.  Energy  capacity 

The  energy  capacity  of  supercapacitor  E22  is  assessed  by 
measuring  the  amount  of  energy  exchanged  (stored  or  released) 
within  the  differential  capacitance  during  the  charging  and  dis¬ 
charging  states.  The  energy-voltage  characteristic  during  25  A 
charging-discharging  is  shown  in  Fig.  11,  together  with  the 
approximation  results  based  on  the  constant  capacitance  method 
and  the  proposed  bivariate  quadratic  method.  In  the  25  A  cycling, 
the  exchanged  energy  is  measured  as  2449.6  J  during  charging  and 
2124.8  J  during  discharging.  In  the  charging  state,  the  estimation  of 
capacitor  energy,  Ec,  and  its  accuracy  with  respect  to  the  recorded 
result  are  2277.8  J  (92.99%)  using  the  conventional  constant 
capacitance  method,  whilst  it  is  2448.2  J  (99.94%)  using  the  author’s 
bivariate  quadratic  method.  In  the  discharging  state,  the  estimation 
of  capacitor  energy  and  its  accuracy  with  respect  to  the  recorded 
result  are  2039.4  J  (95.98%)  using  the  conventional  constant 
capacitance  method,  whilst  it  is  2125.1  J  (100.01%)  using  the  au¬ 
thor’s  bivariate  quadratic  method. 


Fig.  8.  The  charge-voltage  characteristic  of  supercapacitor  E22  under  25  A 
discharging. 
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Fig.  9.  The  differential  capacitance  of  supercapacitor  E22  derived  from  the  bivariate 
quadratic  method  at  different  current  levels. 


Fig.  10.  The  dynamic  voltage  response  of  supercapacitor  E22  under  25  A  cycling. 


Fig.  11.  The  energy-voltage  characteristic  of  supercapacitor  E22  under  25  A  charging¬ 
discharging. 


Table  2 

The  capacitor  energy  and  energy  equivalent  capacitance  of  supercapacitor  E22  at 
different  current  levels. 


Current  /Sc  (A)  Actual  Ec  (J) 

Constant  capacitance 

Bivariate  quadratic 

MJ) 

Ce(  F) 

MJ) 

Ce(F) 

Charging 

25 

2449.6 

2277.8  (92.99%) 

630.41 

2448.2  (99.94%) 

677.58 

50 

2368.0 

2199.4  (92.88%) 

629.21 

2365.7  (99.90%) 

676.80 

100 

2303.7 

2143.2  (93.03%) 

627.20 

2304.6  (100.04%) 

674.43 

200 

2175.4 

2022.6  (92.97%) 

612.44 

2173.0  (99.89%) 

658.00 

Discharging 

25 

2124.8 

2039.4  (95.98%) 

641.53 

2125.1  (100.01%) 

668.49 

50 

2102.1 

2004.1  (95.34%) 

636.23 

2101.4  (99.97%) 

667.12 

100 

1942.7 

1846.4  (95.04%) 

635.79 

1941.9  (99.96%) 

668.70 

200 

1695.1 

1614.9  (95.27%) 

629.57 

1692.0  (99.82%) 

659.64 

Similar  estimation  results  were  obtained  from  constant  current 
cycling  experiments  at  50  A,  100  A  and  200  A,  and  then  compiled 
and  shown  in  Table  2.  In  terms  of  estimation  accuracy,  the  proposed 
bivariate  quadratic  method  achieves  a  relative  error  of  just  0.2% 
compared  to  the  actual  measurements  whilst  the  constant  capac¬ 
itance  method  has  an  error  range  of  4%-7%.  The  energy  equivalent 
capacitance,  Q,  is  the  capacitance  of  a  linear  capacitor  that  contains 
the  same  effective  energy  amount  as  the  differential  capacitance  of 
supercapacitor.  The  CE  is  around  678  F  in  the  charging  state  whilst 
668  F  in  the  discharging  state.  These  values  are  104.2%  and  102.8% 
respectively  of  the  capacitance  value  stated  on  the  nameplate  of 
supercapacitor  E22. 

An  accurate  estimation  of  the  supercapacitor  energy  capacity, 
which  relates  to  the  value  of  the  differential  capacitance,  is 
particularly  important  for  the  estimation  of  the  state  of  charge  of 
supercapacitors.  The  energy  capacity  between  the  half  rated 
voltage  and  the  rated  voltage  is  computed  as  75%  of  the  total  en¬ 
ergy,  based  on  the  constant  capacitance  method.  However,  practi¬ 
cally  this  is  not  the  case  for  supercapacitor  E22  due  to  the  influence 
of  differential  capacitance.  By  referring  to  the  distribution  of  energy 
capacity  plot  in  Fig.  12,  approximately  79%  of  the  device  energy 
capacity  is  available  between  the  half  rated  voltage  and  the  rated 
voltage. 


Fig.  12.  The  distribution  of  energy  capacity  in  supercapacitor  E22  under  25  A  charging. 
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5.  Conclusions 

In  this  paper,  the  authors  have  developed  a  bivariate 
quadratic  approximation  method  and  applied  it  for  the  first  time 
to  represent  the  dynamic  electrical  behaviour  of  supercapacitors 
under  constant  current  cycling  across  an  operating  current 
range  of  25  A-200  A.  This  proposed  method  allows  the  differ¬ 
ential  capacitance  and  energy  capacity  of  supercapacitors  to  be 
derived  directly  from  the  terminal  measurements  of  voltage  and 
current.  This  eliminates  the  need  to  establish  the  relationship  of 
capacitance  versus  voltage  from  multiple  pulsed-current  or 
excitation  signal  tests  at  different  bias  levels  [4,11-14].  The 
approximation  results  of  the  bivariate  quadratic  method  are  in 
close  agreement  with  the  experimental  measurements  at 
different  current  levels  of  25  A-200  A.  For  the  estimation  of  the 
capacitor  energy,  the  relative  error  is  0.2%  using  the  proposed 
method  whilst  in  a  range  of  4%— 7%  using  the  constant  capaci¬ 
tance  method. 

The  quantification  of  differential  capacitance  and  energy  ca¬ 
pacity  leads  to  three  important  observations  in  the  study  of  the 
electrical  characteristics  of  supercapacitors.  Firstly,  the  increase 
of  capacitance  with  respect  to  voltage  is  more  complex  than  a 
linear  function.  Secondly,  the  capacitance  exhibits  variation  be¬ 
tween  charging  and  discharging  states,  and  finally,  approxi¬ 
mately  79%  of  the  energy  capacity  is  available  between  the  half 
rated  voltage  and  the  rated  voltage,  not  75%  as  typically  under¬ 
stood.  The  proposed  method  therefore  can  achieve  a  more  ac¬ 
curate  electrical  model  of  supercapacitors  than  the  conventional 
constant  capacitance  method,  and  it  is  hoped  that  this  new 
technique  will  be  further  adopted  in  the  state  of  charge  estima¬ 
tion  of  supercapacitors,  with  the  aim  to  reduce  the  cumulative 
energy/state  of  charge  error  over  repeated  cycles.  In  practice,  the 
proposed  method  can  be  implemented  as  an  online  diagnosis 
tool  to  analyse  the  health  condition  of  supercapacitors  with  high 
accuracy,  by  using  only  the  terminal  measurements  of  voltage 
and  current. 
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